Ev idence of Convective Redistribution of Carbon Monoxide in AuraTES and MLS Observations 


1 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 


Evidence of Convective Redistribution of Carbon Monoxide in Aura Tropospheric 
Emission Sounder (TES) and Microwave Limb Sounder (MLS) Observations 

M. Many in, A. R. Douglass, M. R. Schoeberl 


Michael Many in 

Science Systems and Applications, Inc., Lanham, Maryland, USA 


Email: Michael.Manvin@nasa.uov 

Address: 

NASA/GSFC 


Laboratory for Atmospheres 
Code 613.3 

Greenbelt, Maryland 20771 


Anne Douglass 

NASA Goddard Space Flight Center, Greenbelt, Maryland, USA 

Email: An ne. R . Doug 1 ass@nasa. gov 

Address: 

NASA/GSFC 

Laboratory for Atmospheres 
Code 613.3 

Greenbelt, Maryland 20771 
Mark Schoeberl 

Science and Technology Corporation, Hampton, Virginia, USA 

Email: Mafk.Schoeberl@mac.com 

Address: 

10005 Old Columbia Rd. 

Suite M- 150 
Columbia, MD 2 1 046 



2 

Evidence of Convective Redistribution of Carbon Monoxide in Aura TES and MtS Observations 

3 1 Abstract 

32 

33 Vertical convective transport is a key element of the tropospheric circulation. 

34 Convection lofts air from the boundary layer into the free troposphere, allowing surface 

35 emissions to travel much further, and altering the rate of chemical processes such as 

36 ozone production. This study uses satellite observations to focus on the convective 

37 transport of CO from the boundary layer to the mid and upper troposphere. Our 

38 hypothesis is that strong convection associated with high rain rate regions leads to a 

39 correlation between mid level and upper level CO amounts. We first test this hypothesis 

40 using the Global Modeling Initiative (GMI) chemistry and transport model. We find the 

4 1 correlation is robust and increases as the precipitation rate (the strength of convection) 

42 increases. We next examine three years of CO profiles from the Tropospheric Emission 

43 Sounder (TES) and Microwave Limb Sounder (MLS) instruments aboard EOS Aura. 

44 Rain rates are taken from the Tropical Rainfall Measuring Mission (TRMM) 3B-42 

45 multi-satellite product. Again we find a correlation between mid-level and upper 

46 tropospheric CO, which increases with rain rate. Our result shows the critical importance 

47 of tropical convection in coupling vertical levels of the troposphere in the transport of 

48 trace gases. The effect is seen most clearly in strong convective regions such as the Inter- 

49 tropical Convergence Zone. 

50 

51 1. Introduction/Background 


52 
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53 Convection plays a vital role in tropospheric circulation as well as the distribution of 

54 trace gases. Surface level trace gas emissions are usually mixed within the atmospheric 

55 boundary layer (BL) via turbulent diffusion but are frequently blocked from entering the 

56 free troposphere (FT) by the inversion layer at the top of the BL. Convection, along with 

57 weather fronts and orographic processes, overcome the BL inversion trapping, moving air 

58 into the FT [Liu et al., 2003], Deep convection reaches the tropopause and sometimes 

59 even overshoots into the lower stratosphere. Cotton et al. [1995J estimated that 

60 convective (cloud related) processes vent the mass of the boundary layer tens of times 

61 over the course of a year. 

62 

63 Once convection lofts constituents into the FT, they can be transported intercontinental 

64 distances because of the consistent large scale flow in that region. Convection also 

65 contributes to inter-hemispheric transport: strong convergence from the Intertropical 

66 Convergence Zone (ITCZ) convective systems can draw air from both hemispheres and 

67 loft it high enough to be pulled into the Hadley circulation. In addition to these effects on 

68 transport, convection can also dramatically alter the rate of chemical processes because of 

69 its speed. Convective time scales are very short, on the order of hours for surface to UT, 

70 as opposed to days for synoptic scale processes. The rapid transport of chemical species 

71 means that short-lived species can still reach the upper troposphere. For example, when 

72 urban air with NO x and other ozone precursors is quickly lofted to the FT, ozone 

73 production is greatly enhanced [Pickering, 1992] because the transport is fast compared 

74 to the NO x lifetime. Convective transport also extends the lifetime of some species by 
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75 reducing their exposure to loss processes in the lower troposphere; this expands their 

76 range of chemical influence [ Dickerson et al., 1987]. 

77 

78 Convective processes are complex, involving surface radiative heating, updrafts from 

79 lower tropospheric levels, detrainment at heights up to and above the tropopause, 

80 condensation, latent heat release, precipitation, downdrafts and evaporative cooling (e.g. 

81 Gamache [1982]). Other factors include melting ice particles, and turbulent mixing due to 

82 vertical shear and gravity wave instabilities. The scale of convective systems can range 

83 from isolated thunderstorms up to mesoscale complexes measuring 1000+ km. The 

84 wide-ranging impact of convection on chemical transport makes it crucial for us to 

85 quantify convective transport and evaluate its representation in models. 

86 

87 Previous studies of convection have focused on various features, including (i) showing 

88 evidence of convective transport from BL to upper troposphere/lower stratosphere region 

89 [Poulida, 1996; Randel and Park, 2006; Li et al., 2005; Park et al., 2007; Ricaud et al., 

90 2007; Ziemke et al., 2009] (ii) determining the height of convective outflow [Gettelman 

91 and Forester, 2002; Froyd et al., 2009] (iii) identifying cases where outflow penetrates 

92 the tropopause, into the LS [Hegglin, 2004; Ray et al., 2004; Liu and Zipser, 2005; 

93 Chaboureau et al., 2007] (iv) comparing post-convective UT values measured by 

94 satellite to model predictions [Liu et al., 2003; Holland et al., 2009; Mari et al., 2000] 

95 and (v) the diurnal cycle of convection [Tian et al., 2004; Liu and Zipser, 2009], 


96 
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97 One result which is especially pertinent to this study is the trimodal distribution of 

98 tropical clouds, which can be traced to levels of convective detrainment. Johnson el al. 

99 [1999] reported clouds forming at three main levels, as detected by shipbome radar 

100 during the TOGA COARE experiment: low cumulus near 2 km and cumulonimbus at 1 5 

101 to 16 km were most prominent, but cumulus congestus at approximately 5 km were 

102 frequent as well. These levels correspond well with stable layers associated with the trade 

103 wind inversion, the tropopause, and the 0°C melting level isotherm. Such stable levels are 

104 understood to be preferential regions of convective outflow (e.g. Zuidema [1998]). The 

105 mid-level cloud layer has also been documented by Dessler et al. [2006] using GLAS 

106 data, and by Haynes and Stephens [2007] and Riley and Mapes [2009] with CloudSat 

107 data. We will examine trace gas concentrations at the middle and upper levels of the 

1 08 trimodal distribution, to demonstrate the effects of convective transport. 

109 

1 10 The goal of this study is to use satellite measurements of CO combined with satellite- 

1 1 1 based locations of convective systems to quantify the impact of the convective systems 

112 on trace gas distributions. This study uses a 3 year archive of satellite data to search for 

1 1 3 evidence of convective transport of CO from the BL to mid and upper tropospheric 

1 14 levels. Observations are restricted to the tropics, where incoming solar radiation is 

1 15 strongest and convection is frequent. We use CO retrievals from the Tropospheric 

1 16 Emission Sounder (TES) on EOS Aura to detect mid-level concentrations, and CO 

1 1 7 retrievals from the Microwave Limb Sounder (MLS) aboard Aura for upper tropospheric 

118 values. To identify areas of convective transport we use estimates of rain rate derived 
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1 19 from microwave and IR satellite measurements, calibrated with Tropical Rainfall 

120 Measuring Mission (TRMM) observations, as a proxy for convective strength. 

121 

122 Carbon monoxide is a product of fossil fuel combustion and biomass burning, as well as 

123 the oxidation of hydrocarbons, particularly methane [Jacob, 1999], As such, its highest 

124 concentrations are near the Earth’s surface. CO is lost in the troposphere primarily by 

125 reaction with the hydroxyl radical OH. The lifetime of CO generally ranges between one 

126 and two months, depending on season and location within the atmosphere. Carbon 

127 monoxide is not soluble and therefore not subject to rainout. These qualities make it an 

128 excellent tracer. CO is monitored for many reasons, among which: it is a significant 

129 anthropogenic pollutant, it is a precursor to O 3 in the troposphere, and it is a sink for OH 

130 (thus modulating the oxidation of many greenhouse gases). Global CO measurements are 

131 currently provided by a number of space-borne instruments, including Measurements Of 

132 Pollution In The Troposphere (MOPITT) [Drummond, 1992] aboard the EOS Terra 

133 satellite, Atmospheric Infrared Sounder (AIRS) [McMillan et al., 2005] aboard EOS 

134 Aqua, MLS [Waters et al., 2006] and TES [Beer et al., 2001], both on EOS Aura. We 

135 chose to use TES and MLS in this study: TES for its sensitivity in the mid-troposphere 

136 and its small footprint, allowing for unobscured profiles near or between clouds; MLS for 

137 its sensitivity near the tropopause, its capability to see through clouds, and its consistent 

138 proximity to TES profiles (as compared with MOPITT). 

139 

140 To locate convection, one must detect an accompanying phenomenon, such as clouds 

141 resulting from condensing of lofted water vapor, or rainfall from those clouds. Numerous 
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142 approaches to cloud characterization exist, including diagnosing outgoing longwave 

143 radiation (OLR) [ Liebmann and Smith, 1996; Lee et al., 2007], GOES Precipitation Index 

144 (GPI) \Arkin and Meisner, 1987] based on microwave brightness temperatures, and 

145 International Satellite Cloud Climatology Project (ISCCP) cloud classification [Rossow 

146 and Schiffer, 1991] based on cloud top pressure and optical depth. Techniques such as 

147 OLR and ISCCP yield high resolution products, but global coverage is only once or twice 

148 a day. GPI and its derivatives use data from geosynchronous satellites for good spatial 

149 and temporal resolution, but features such as high cirrus and non-raining portions of 

150 cloud anvils pose difficulties. 

151 

152 Microwave instruments that can detect precipitation beneath clouds are quite effective at 

153 locating convective rainfall. Successful estimation of rain rate from passive microwave 

1 54 satellite data was first demonstrated with ESMR on Nimbus satellites [ Wilheit et al., 

155 1973], Current sensors include SSM/I on DMSP spacecraft, and TMI and PR aboard 

156 TRMM. Product resolution is quite good, but coverage is incomplete. 

157 Algorithms such as the TRMM 3B-42 (see section 2.2) estimate rain rates with 

158 microwave data where available and calibrated IR data elsewhere, in addition to rain 

159 gauge values. This provides full coverage in the low and mid latitudes, every 3 hours, at 

160 quarter degree resolution. We presume that higher rain rates from the TRMM 3B-42 

161 product imply greater likelihood of convection. 

162 

163 2. Approach 


1 64 



Evidence of Convective Redistribution of Carbon Monoxide in AuraTES and MLS Observations 


8 


1 65 The coupling of the BL and FT by convection should have a definite signature in the 

1 66 trace gases. CO should appear aloft following a convective event that ventilates the BL, 

167 and outflow levels in the mid and upper troposphere should both show evidence of this 

168 lofting transport, since detrainment preferentially occurs near the stable layers of 0°C 

169 (occurring between 500 and 600 hPa) and the tropopause (-200 hPa). Thus there should 

1 70 be higher correlation between trace gas concentrations at mid and upper tropospheric 

1 7 1 levels following periods of local convective activity. 

172 

1 73 Our approach is to examine correlations between concentrations of CO in the middle and 

174 upper troposphere, and relate them to proximity to significant rainfall. First we examine 

175 this relationship within the context of a numerical model, the Global Modeling Initiative 

176 (GMI) chemistry and transport model (CTM). Using this model will allow us to examine 

1 77 what might be observed at the model scale of 2.5° longitude x 2° latitude. Then we 

178 search for a similar relationship using data from TES and MLS on Aura. We use TES 

179 pressure level 562 hPa and MLS level 215 hPa, because these correspond well with 

1 80 convective detrainment levels put forth in the trimodal cloud characterization [Johnson et 

181 al., 1999], We use the TRMM-calibrated 3B-42 rain rate product to identify likely 

1 82 instances of convection. 

183 

1 84 One challenge of using TES observations is that moist convection is accompanied by 

1 85 clouds, which prevent infrared-based remote sensing of the lower atmosphere. We 

186 anticipated however that cloud-free observations of the middle troposphere would be 

1 87 available immediately adjacent to convective areas, e.g. between clouds or along the 
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1 88 edges of convective storms. The small TES footprint makes this approach reasonable. 

1 89 We anticipate only a small fraction of the -650 measurements made in the tropics (1 8°N- 

1 90 S) during each global survey will meet our criteria for cloud amount and proximity to 

191 convection. We compensate for infrequency by drawing upon a multi-year archive. 

192 

193 The remainder of the paper is structured as follows: The models (CTM and AGCM) are 

1 94 discussed in section 2. 1 ; the data sources for CO and rain rate are presented in 2.2. 

195 Results are given in section 3, and conclusions follow in section 4. 

196 

197 2.1 GM1 and GEOS-4 

198 

199 NASA’s Global Modeling Initiative has developed a ‘Combo’ CTM that combines 

200 photochemical schemes that are appropriate for the stratosphere and troposphere. The 

201 Combo CTM may be driven by meteorological fields from a general circulation model 

202 (GCM) or from a data assimilation system (DAS). GMI Combo includes gas phase 

203 reactions, photolysis, wet scavenging, dry deposition, heterogeneous chemistry (PSCs) 

204 and emissions [Strahan et ctL, 2007; Duncan et al., 2007], For this study we have used 

205 the “Aura4” GMI experiment, run at 2° latitude by 2.5° longitude, with 55 vertical levels 

206 (hybrid sigma-pressure) up to 0.015 hPa. We use simulated CO that is recorded at times 

207 corresponding to the Aura afternoon overpass times, between 1 and 2 p.m. local time to 

208 best match the observation times of TES on Aura. 


209 
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210 For the Aura4 experiment, GM1 ‘Combo' used meteorological fields taken from the 

2 1 1 Goddard Earth Observing System Data Assimilation System (GEOS DAS) version 4 

212 [Bloom et al., 2005], The GEOS DAS fields are produced at 1° latitude by 1 .25° 

213 longitude on 55 levels; the fields are adapted for use at 2° by 2.5° in the CTM. The 

214 physics parameterization scheme consists of parameterizations for precipitation, 

215 clouds/radiation, turbulent mixing, and surface processes. The deep moist convection 

2 1 6 scheme is that of Zhang and McFarlane [1995], which uses a plume ensemble approach 

217 with convective updrafts and related downdrafts. Shallow convection is covered by the 

218 Hack [ 1 994] scheme. 

219 

220 GEOS-4 DAS assimilates observations from many sources, including SSM/I (total 

221 precipitable water), MODIS (winds), geosynchronous satellites (winds), TOVS sounders 

222 (temperature and moisture), ERS-2 and QuikScat (sea-level winds), along with 

223 rawinsondes (temperature, winds, moisture), and surface stations (temperature, pressure, 

224 winds). The DAS outputs dozens of fields, including those that GMI requires, such as 

225 winds, temperature, precipitation and convective mass flux. Fields are saved at three 

226 hour intervals, in instantaneous and time-averaged forms; GMI ‘Combo’ uses the latter, 

227 horizontally averaged to 2° x 2.5°. In order to identify grid boxes with convective 

228 activity, we used the GEOS-4 DAS output field “PRECON”, which stores the convective 

229 precipitation rate. Each GEOS-4 output file covers a single synoptic time period. In 

230 order to match the local ascending Aura overpass times implicit in each GMI output file, 

23 1 the convective precipitation field had to be extracted from multiple GEOS-4 files 

232 covering 24 hours of simulation. 
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233 

234 2.2 TES, MLS and TRMM (3B-42) 

235 

236 The Tropospheric Emission Spectrometer (TES) is a Fourier transform spectrometer 

237 [ Beer et al., 2001; Beer, 2006]. It is one of several instruments aboard NASA’s Aura 

238 satellite, positioned in an ascending sun synchronous orbit with equator crossing at about 

239 13:45 local solar time. TES measures in the infrared to detect trace gases including CO, 

240 water vapor, O 3 and methane. This study uses nadir profiles acquired during TES Global 

241 Surveys. Nadir view operates with a 5.3 x 8.3 km footprint, and each profile in a Level 2 

242 file consists of data on 65 pressure levels ranging from nominal surface pressure to 0. 1 

243 hPa. Compared to the GMI model, TES has a much higher horizontal resolution, but of 

244 course the TES observations are much sparser than the GMI output (e.g. 9800 tropical 

245 profiles in July 2007, compared with 84,8 1 6 GMI gridbox values). And in terms of 

246 vertical resolution, the model levels within a grid column are independent compared to 

247 low vertical resolution retrieved TES profiles. 

248 

249 In the tropics, a TES tropospheric CO profile typically has between 1 and 2 degrees of 

250 freedom (DOF), a measure of the number of “independent” pieces of information that 

25 1 contributed to that profile. These pieces of information are used to modify the a priori 

252 (first guess) trace gas profile. (See Beer [2006] for an overview.) With a greater number 

253 of DOFs it is possible to achieve more variation from the a priori. In December 2005, 

254 TES implemented an Optical Bench warmup which resulted in better instrument 
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255 alignment, and average DOF values for TES CO measurements were significantly 

256 improved. For this reason, we have only used data from January 2006 forward. 

257 

258 Asa consequence of variations in TES DOF, some measurements reflect the actual state 

259 of the atmosphere, while others rely heavily on the a priori. To distinguish between such 

260 cases we consult the averaging kernel diagonal (AKD) values; where these are high 

261 (close to 0.1), there is greater confidence that the true atmospheric values are retrieved 

262 [Halland et al., 2009], We tested several AKD thresholds and use an AKD threshold of 

263 0.085 for most of this paper. In the final figure we show results using two other 

264 thresholds for the AKD. 

265 

266 TES profiles used in this study were Tevel 2, Version 4. They were filtered using the 

267 quality control fields, as in Table 5-3 of the EOS TES L2 Data User ’s Guide [ Osterman 

268 et al., 2008]. In addition, since the TES IR-based retrievals are subject to interference by 

269 clouds, profiles with nominal cloud top pressures less than 572 hPa were discarded. We 

270 used profiles between 1 8°N and 1 8°S because DOFs changed abruptly beyond those 

271 boundaries due to changes in the a priori. 

272 

273 The MLS instrument aboard Aura measures profiles of radiance from the Earth’s limb 

274 from upper troposphere to the mesosphere [Waters, 2006], MLS scans columns ahead of 

275 the satellite, every 1 65km along its orbit. MLS measures CO [Livesey et al., 2006; 

276 Pumphrey et al., 2007], along with a dozen other atmospheric trace gases. As a 
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277 microwave instrument, it can detect constituents despite the presence of aerosols, cirrus 

278 or polar stratospheric clouds. 

279 

280 In this study we use CO as reported by MLS (version 2.2) at 2 1 5 hPa. At this level, the 

28 1 CO morphology and variability are adequately measured, but the values are biased by 

282 approximately a factor of 2 [Livesey et al., 2008], We compensate by dividing all values 

283 in half; we also processed the data according to guidelines in the MLS data quality 

284 document [Livesey et al., 2007]. Each MLS CO profile was paired with the closest TES 

285 profile; the average separation was 195 km, and the average time difference was 

286 approximately 7 minutes. 

287 

288 For rain rate data, we turned to the 3B-42 product from the Tropical Rainfall Measuring 

289 Mission (TRMM) [Simpson et al., 1996]. This mission, a joint effort of the United States 

290 and Japan, was launched in 1997. It supports several instruments including the TRMM 

291 Microwave Imager (TMI), the Precipitation Radar (PR), and the Visible and Infrared 

292 Scanner (VIRS) [Kummerow et al., 1998]. The 3B-42 algorithm [Huffman et al., 2007] 

293 combines observations from several satellites, including TRMM, with rain gauge 

294 analyses. Computation of the multi-satellite portion begins with a merged microwave 

295 product including TMI and SSM/I rain rate estimates. Then VIRS data is used to produce 

296 an infrared/rain rate relationship, which is applied to geosynchronous IR data, such as 

297 that from the Geostationary Operational Environmental Satellite (GOES) series. The 

298 microwave products are high resolution (e.g. TMI is 5.1 km), but measurements 

299 accumulated over timescales of up to three hours do not completely cover the tropics. 
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300 The IR measurements are similar resolution (e.g. GOES is 4 km), and a merged product 

301 from several geo-lR satellites provides excellent coverage (60° N-S) at half hour intervals 

302 [Janowiak et al. , 2001]. Incorporating both microwave and IR data, the 3B-42 rain rate 

303 product is generated for latitudes covered by TRMM (50° N-S) at quarter degree 

304 resolution and 3 hour intervals. 

305 

306 3B-42 does not distinguish between convective rain, which is linked to strong vertical 

307 transport, and stratiform rain, which is not. However, it is generally true that greater rain 

308 rates are more likely associated with convective systems. Indeed, the ratio of convective 

309 to stratiform rain rates in the tropics is 4.1, when measured at PR resolution [Schumacher 

310 andHouze, 2003]. So the error introduced by the inclusion of stratiform rain rates in this 

3 1 1 study should be relatively small. 

312 

313 Rain rates were used as follows: For a given pair of TES and MLS CO profiles we 

314 searched a neighborhood of 3B-42 grid boxes centered on the spot midway between the 

315 profiles. Since 3B-42 rain rates are available every three hours, we used the values from 

316 the 3-hour window prior to the TES observation, to ensure that we would see the effects 

3 1 7 of vertical mixing. 

318 

319 3. Results 

320 

321 3.1 GMI and GEOS-4 


322 
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323 We first analyze simulated upper and mid-tropospheric CO to demonstrate the expected 

324 relationship between composition and convective transport as implemented in the GMI 

325 CTM. 

326 

327 3.1.1 Scatter plots varying rain rate threshold 

328 

329 Daily GMI CO concentrations from July 2007 were compared at two different pressure 

330 levels, 562 hPa and 215 hPa, selected to match the TES and MLS pressure levels used in 

331 the next section. The 562 hPa level represents the level of expected convective 

332 detrainment around the 0°C isotherm, while 2 1 5 hPa corresponds with outflow near the 

333 tropopause. As shown in the first scatter plot of Figure 1 , the majority of CO values on 

334 both pressure levels range between 50 and 1 50 ppbv. Panel (a) includes values from all 

335 GMI grid boxes in the tropics, 19°S - 19°N, regardless of rain rate. The broad pattern 

336 matches our expectations: Points along the diagonal are consistent with vertical 

337 convective transport to both levels of the troposphere; areas with substantial CO at lower 

338 altitudes but minimal amounts above may represent cases where shallow convection has 

339 deposited CO emissions from the boundary layer up into the mid-troposphere. Areas 

340 with greater CO concentrations at 21 5 hPa may represent cases where deep convection 

341 has been dominant. The off-diagonal cases could also be due to the effects of horizontal 

342 advection transporting CO from other convective areas. 

343 

344 Figure 1 , panels (b)-(d) show subsets of the original data, where profiles were discarded 

345 if the convective rain rate was less than 0.5, 1 .0 and 2.0 mm/hr respectively. As 
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346 convective activity increases from one panel to the next, profiles with markedly different 

347 CO concentrations at the two levels are eliminated. This fits with our hypothesis, in that 

348 convective transport is coupling the middle and upper levels of the troposphere. In the 

349 last panel, where convective rain rates are highest, the remaining points are clustered 

350 along the diagonal, and the correlation coefficient is quite high, at 0.967. This simple test 

351 demonstrates that the GMI implementation of trace gas transport, as driven by GEOS-4 

352 winds and convective mass flux, is consistent with our conceptual model of the effect 

353 vertical convective mixing will have on the trace gas distributions. This is an idealized 

354 case with which to compare the TES/MLS observations. 

355 

356 3.1.2 Geographic Location of GEOS-4 Rain 

357 

358 Another useful plot for comparing the model and satellite results is the geographic 

359 location of dominant rainfall. The areas of convective rain in the GEOS-4 DAS are 

360 shown in Figure 2. The grid boxes are shaded based on the number of days that met or 

361 exceeded the rain rate threshold. Note the prominence of the ITCZ in panels (a) and (b), 

362 and the areas of most frequent convective rain such as Costa Rica and north of Papua 

363 New Guinea, in panel (c). These plots confirm the success of GEOS-4 DAS in 

364 reproducing realistic rain distribution patterns. 

365 

366 3.2 TES, MLS and TRMM (3B-42) 


367 
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368 The TES and MLS observations are quite sparse in comparison to model output. Within 

369 our constraints of tropical latitudes and sufficiently cloud-free observations, a month of 

370 TES Global Survey data yields roughly 30 to 60 profiles. Filtering these using rain rate 

37 1 thresholds quickly reduces the sample size to numbers too small for a meaningful 

372 analysis. So we chose to include three full years of TES and MLS data in this study, 

373 from 2006 through 2008. 

374 

375 In order to estimate the strength of convection in the neighborhood of TES and MLS 

376 measurements, we consulted the TRMM-calibrated rain rate data for the 3-hour interval 

377 prior to the TES overpass time. Our initial tests included an area similar in size to a GMI 

378 gridbox: we searched an 11x7 patch of TRMM values centered between the TES and 

379 MLS profiles, and took the maximum rain rate. (Since TRMM 3B-42 is 0.25° resolution, 

380 the area was 2.75° in longitude by 1.75° in latitude.) In addition we tested our approach 

381 with smaller patches of TRMM grid boxes, including 7x5, 5x3 and 3x3. As the area in 

382 which we searched for rain grew smaller, the effects of the vertical mixing became more 

383 apparent. So we begin by examining the 3x3 case, wherein each patch is a square of 9 

384 grid boxes; the maximum rain rate is sought within a 0.75° x 0.75° region. 

385 

386 3.2. 1 Scatter plots varying rain rate thresholds 

387 

388 Plots of CO correlation between pressure levels 562 and 215 hPa are shown in Figure 3. 

389 Panel (a) includes all qualifying profiles, regardless of local rain rate. It shows that the 

390 majority of values at midlevel are between 50 and 150 ppbv, while the range on the upper 
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39 1 level is between 50 and 125 ppbv. The (solid) line of linear fit shows that greater CO 

392 concentrations exist at lower altitudes. This can be attributed to diffusive transport from 

393 the boundary layer, particularly during times of weak trade inversion. It may also be due 

394 in part to the approximate correction we have applied to compensate for the MLS high 

395 bias. Greater CO concentrations at mid-level were also present in GMI (see the line of 

396 linear fit in Figure 1 ), but in the TES/MLS case this is much more pronounced. 

397 

398 Panels (b)-(d) in Figure 3 show the relationship between middle and upper tropospheric 

399 CO as the rain rate requirement for the TES/MLS observation neighborhood is increased. 

400 As the threshold is increased from 1 up to 3 mm/hr, the correlation steadily increases as 

401 well. TES/MLS values associated with higher rain rates are clustered around the 

402 diagonal, as was the case in the GMI plots (Figure 1). Clustering is tighter for GMI, 

403 because values are reported immediately after the convective mass flux is applied to the 

404 trace gas field, as opposed to the TES/MLS measurements which were taken up to 6 

405 hours after rainfall. The interval between convective rain and CO observations gives 

406 advection time to alter distributions. 

407 

408 Part of the difficulty in characterizing patterns in this data is the remarkably low number 

409 of profiles satisfying our criteria, as shown in Table 1 . This is due in part to the rejection 

410 of cloudy TES profiles. For the 0 mm/hr case, 84% of the profiles were filtered out based 

411 on their cloud top pressure. In the case of 3 mm/hr, 90% were removed. The linkage of 

412 additional clouds with higher rain rate was expected, since convective rain is 

413 accompanied by updrafts, leading to detrainment, condensation, and therefore clouds. 
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414 

415 To establish a level of confidence in the correlation values computed for small sample 

416 sizes, we employed a bootstrap approach: we computed the correlation for a randomly 

417 chosen subset of half the points, repeated this process 1000 times, and made histograms 

418 of the results. We accepted the correlation over the full sample only when it matched the 

419 average of the bootstrap correlations to within 0.05, and when the central 90% of the 

420 bootstrap values did not exceed a range of 0.5. Each of the correlations listed in Table 1 

421 satisfied these criteria. The bootstrap histograms are shown as insets in Figure 3, with 

422 markers denoting the 5%, 95% and average values. 

423 

424 3.2.2 Geographic Location of TES/MLS Observations 

425 

426 Because strong convection is usually accompanied by substantial clouds, we anticipated 

427 that non-cloudy TES profiles linked with higher rain rates would be adjacent to regions of 

428 heavy rain, rather than within such regions. The locations of accepted TES/MLS profile 

429 pairs are shown in Figure 4. In panel (a), rain rates are not considered, and the pattern 

430 arises solely from the elimination of cloudy profiles. The areas of highest number of 

43 1 observations are those that tend to have low levels of rain, such as Saudi Arabia and the 

432 horn of Africa. The fewest profiles are found in areas of high rainfall, such as the ITCZ, 

433 Maritime Continent, Indian Ocean and Brazil. So, as expected, most non-cloudy TES 

434 observations are in the areas of least convective rain. Once a rain rate threshold is 

435 applied a different pattern emerges. (Note that the color scale for panels (b) and (c) is 

436 different from that in the first panel.) As anticipated, the profiles in rainy neighborhoods 
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437 tend to be located adjacent to areas of frequent rain, such as the ITCZ, Brazil, and the 

438 Indian Ocean. We interpret these to be instances where TES measures CO in sufficiently 

439 cloud-free columns along the edges of convective storm systems. The convective lofting 

440 of CO by such systems is evident in the TES and MLS data. With higher rain rates there 

441 is greater likelihood of convective activity in the area, and therefore greater vertical 

442 mixing of CO and higher correlation between levels. 

443 

444 3.2.3 Correlation affected by size of TRMM neighborhood, and AK Diagonal 

445 

446 How does the correlation between mid and upper levels of CO change as the distance 

447 from convection increases? And how is the correlation affected by our choice of AKD 

448 threshold? To answer these questions, we performed two sensitivity experiments. First 

449 we varied the size of the search neighborhood in the TRMM-based rain rate data, testing 

450 four different patch sizes: 1 1x7, 7x5, 5x3 and 3x3 grid boxes, in the 0.25° x 0.25° grid. 

45 1 Larger search areas allow for rain rate maxima to be identified at greater distances from 

452 the CO profiles. In the second sensitivity test we varied the minimum TES AKD 

453 accepted: 0.08, 0.085 and 0.09. This was done to determine the trend in correlations as 

454 TES values gained greater independence from the a priori and more accurately reflected 

455 the state of the atmosphere. The results from both of these experiments are shown in 

456 Figure 5. Multiple patch sizes are plotted in each panel, and a different minimum AKD 

457 applies to each panel. (Note that the original series of correlations, recorded in Table 1 , 

458 appears as the sequence of red asterisks in the middle panel.) 


459 
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460 The influence of TRMM search area size is evident in panel (a). Each gray square 

461 designates the CO correlation derived from TES/MLS profile pairs close to a sufficiently 

462 high TRMM rain rate, determined by searching an 11x7 neighborhood of TRMM values 

463 (2.75° x 1 .25°). The blue triangles show the results from searching smaller 

464 neighborhoods of 7x5 grid boxes (1 .75° x 1 .25°). The green diamonds reflect 5x3 patches 

465 (1.25° x 0.75°), and the red asterisks use 3x3 patches (0.75° x 0.75°). 

466 As the search area is constricted, the average distance between the TES/MLS profiles and 

467 maximum rain rate is reduced, so the lofting effects of convection should be more evident 

468 in the CO measurements, and we anticipated higher correlation as a result. This is clearly 

469 exhibited in the plots; while the 1 1x7 sequence does not show an increase in correlation 

470 as local maximum rain rate increases, the 7x5 sequence does show this relationship, and 

471 the 5x3 and 3x3 cases even more so. This demonstrates that trace gas correlation between 

472 the middle and upper troposphere depends on proximity to convection as well as 

473 convective intensity. 

474 

475 The series of three panels shows the impact of more stringent TES AKD requirements. 

476 Two effects are apparent. First, as average AKD increases, the CO correlation for a given 

477 rain rate tends to increase as well. This means that the trends we are seeing derive 

478 primarily from the atmospheric state, and should not be attributed to the TES a priori. 

479 The second effect is the reduction in qualified TES profiles, resulting in smaller sample 

480 sizes over which to evaluate correlations. This paucity of data is further exacerbated 

48 1 when profiles are filtered based on rain rate. With TES AKD > 0.08, before rain rate is 

482 considered, there are 17,358 profile pairs from the 2006-2008 period; with AKD > 0.085 
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483 there are 9,945 pairs, and with AKD > 0.09 only 5,091 pairs qualify. Filtering these sets 

484 with increasing rain rate thresholds quickly reduces the sample sizes to fewer than 50, 

485 and soon thereafter the bootstrap confidence test (see section 3.2.1) begins to fail, 

486 discontinuing each sequence. For example, the 3x3 sequence in panel (c) is computed for 

487 rain rates at or above 0, 0.5 and 1.0 mm/hr, using 5,091, 83 and 41 profile pairs 

488 respectively. At 1.5 mm/hr, too few profiles remain to calculate a statistically meaningful 

489 correlation. We assert that, despite the short series of correlations plotted in the last panel, 

490 they confirm the trends seen in the prior two panels, namely that the proximity and 

491 intensity of convection influences the degree of correlation in middle and upper 

492 tropospheric concentrations of CO. 

493 

494 4. Conclusion 

495 

496 Using rain activity as a proxy for convection we have demonstrated a link between 

497 tropical convection and CO concentrations in the middle and upper troposphere, at levels 

498 corresponding to the stable layers at 0°C and the tropopause. In other words, convective 

499 activity is pumping CO from the surface into the mid and upper troposphere - chemically 

500 coupling the tropospheric layers. This association was first demonstrated in the GMI 

501 CTM. Correlations between CO concentrations at 562 and 215 hPa increased as the 

502 convective precipitation rate increased. This conformed to our expectations. In addition, 

503 the locations of high convective rain rates for the month we studied aligned well with 

504 standard patterns of tropical rain, e.g. the ITCZ. 

505 
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506 For the observational study, CO measurements from TES at 562 hPa (unobscured by 

507 clouds) and from MLS at 215 hPa were examined, together with rain rates from the 

508 TRMM 3B-42 product. First we examined the influence of rain rate on the correlation, 

509 using a 0.75 x 0.75 search area within TRMM data for each pair of TES/MLS profiles, 

510 and using TES profiles with at least an AKD value of 0.085. As the rain rate threshold 

511 was increased, the correlation between CO burdens at the two levels grew stronger, as we 

512 anticipated; greater convective vertical mixing is concomitant with heavier rain. 

513 Flowever, the sample size (profile count) fell dramatically as we restricted ourselves to 

514 higher local rain rates. This was because rain is accompanied by clouds (convective rain 

515 leads to additional cloud formation), and we discarded TES profiles that were not 

5 1 6 sufficiently cloud-free. 

517 

5 1 8 The influence of clouds was also evident in the location of qualifying profiles. Most 

5 1 9 profiles clustered in areas of minimal rain, such as the horn of Africa. When rain rate 

520 thresholds were enforced, the remaining profiles appeared at the periphery of major 

521 cloudy areas such as the ITCZ and Brazil. These plots differed greatly from the GMI 

522 location plots, principally due to the cloud-free condition we imposed on TES profiles. 

523 Fortunately, the small footprint of TES permitted some measurements not obscured by 

524 clouds immediately adjacent to regions with rain. We extended the study to include 36 

525 months of data to provide a sufficient number of non-cloudy profiles. Also we used a 

526 bootstrap statistical method to gauge confidence in cases with only a small number of CO 

527 measurements. By these means we were able to show linkage between rain rate and CO 

528 burden correlation similar to that shown with the model. 
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529 

530 We explored the sensitivity of this result to two parameters: the rain search area, and the 

531 TES averaging kernel diagonal values. For the former, we varied the extent of search 

532 within TR.MM 3B-42 for a maximum rain rate to associate with each TES/ML.S profile 

533 pair. TRMM data neighborhoods of four sizes, ranging from 2.75°xl .75° (approximately 

534 the size of one GMI gridbox), to 0.75°x0.75° were tested. For each areal size, our first 

535 result was reproduced: increased rain rates led to higher correlations. Smaller search 

536 domains also tended to yield better correlations, particularly for greater rain rates. This 

537 result is consistent with the idea that closer proximity to convective mixing translates into 

538 a stronger observable signal. When we required higher AKD values for TES, the 

539 correlations again tended to rise, giving us confidence that the signal in CO reflected 

540 actual atmospheric concentrations, rather than strong influence from the a priori. We 

541 showed significant correlations up to AKD threshold 0.09, beyond which the profile 

542 count was too small for meaningful results. 

543 

544 We have demonstrated that tropical convection plays an important role in transporting 

545 trace gases vertically within the troposphere, with effects in keeping with the trimodal 

546 characterization of tropical clouds, in which boundary layer parcels are convectively 

547 lofted to stable levels in the troposphere. Satellite observations show that this process 

548 operates on a wide scale throughout the tropical regions where strong convection is 

549 prevalent. 


550 
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551 In this paper we examined how correlation was affected by proximity to rainfall but 

552 subsequent advection of the CO plume was not considered. A follow-up study could 

553 focus on cases where convection was located upwind of the CO measurement. 

554 

555 
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778 Figure 1 . GM1 CO concentrations at middle and upper tropospheric levels, using one 

779 month of daily output between 19°S-19°N. Panel (a) shows values from all grid boxes. 

780 In subsequent panels, values from a grid box are only included if the convective 

78 1 precipitation rate within that grid box meets the given threshold. The correlation of CO 

782 between the two levels increases as the rain rate increases. 

783 

784 Figure 2. Each grid box shows the number of days (day count) in July 2007 for which 

785 GEOS-4 computed a tropical convective rain rate satisfying the given threshold, at Aura 

786 afternoon overpass times. Areas of most frequent convective rain correspond with well- 

787 known areas of convection e.g. the ITCZ and the Maritime Continent. In panel (a), rain 

788 rate > 0.05 mm/hr, panel (b) > 0.5 mm/hr, and panel (c) > 1 mm/hr. 

789 

790 Figure 3. CO concentrations at middle vs. upper tropospheric levels, as measured in the 

791 tropics (18°S-18°N) from Jan 2006 through Dec 2008. TES Global Survey data at 562 

792 hPa were used, only where the Averaging Kernel diagonal was at least 0.085. MLS data 

793 were used at 2 1 5 hPa, with values halved to compensate for a factor of 2 bias at this level. 

794 The profiles were screened for quality and cloud presence as described in the text. Panel 

795 (a) shows values from all qualifying profiles. In subsequent panels, a profile is only 

796 included if the maximum local TRMM-based rain rate (within a 3x3 grid patch) meets the 

797 given threshold. As the 3B-42 rain rate threshold increases, so too does the con-elation of 

798 CO between the two levels; at the same time, the number of points decreases pretty- 

799 dramatically. Histogram insets show the distribution of 1000 bootstrap test correlations. 
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800 ranging from -1 to 1 . Markers below the histograms show the 5%, average, and 95% 

801 locations; markers above show the correlation coefficients using all data. 

802 

803 Figure 4. Similar to Figure 2; however, this figure shows the geographic distribution of 

804 TES/MLS profile pairs with TRMM rain rates satisfying thresholds 0.0, 0.05 and 1.0 

805 mm/hr (panels (a) and (c) correspond to panels (a) and (b) in Figure 3). The profiles here 

806 are binned at 2° x 2.5° to match Figure 2. In panel (a), no rain rate threshold is applied, 

807 and the cloud-limiting filter for TES leaves clusters of profiles in areas of very little rain 

808 (e.g. parts of Africa, Saudi Arabia and the South Atlantic). Panels (b) and (c) show the 

809 effects of imposing rain rate thresholds of 0.05 and 1 .0 mm/hr; this greatly reduces the 

810 number of profiles, so a different color scale is used. As the required rain rate increases, 

81 1 the remaining profiles are frequently located near areas of prominent rain (e.g. the ITCZ, 

812 Maritime Continent, and Brazil). 

813 

814 Figure 5. Correlation of CO measurements from TES (562 hPa) and MLS (215 hPa), 

815 plotted against rain rate thresholds satisfied by TRMM 3B-42. Four methods of rain rate 

816 determination are shown in each panel; they differ only in the number of TRMM values 

817 consulted in the vicinity of the TES/MLS observations, ranging from 77 (1 1x7) down to 

818 9 (3x3). Each panel imposes a different AKD minimum for the TES data: 0.08, 0.085 and 

819 0.09. Data spans three full years, 2006-2008. The points on the red/asterisk line (3x3 

820 patches) in panel (b) correspond to the correlation coefficients shown in Table 1. 

821 Correlation is not recorded if it does not satisfy the bootstrap criteria described in section 

822 3.2.1. 
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823 

824 Tables 

825 

826 


Rain Rate 
Threshold 
(mm/hr) 

0 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

Correlation 

Coefficient 

0.454 

0.538 

0.604 

0.618 

0.719 

0.780 

0.804 

Profile count 

9940 

213 

108 

64 

45 

33 

22 


827 

828 Table 1. Interlevel CO Correlation Coefficients 

829 CO correlations between middle and upper tropospheric levels associated with various 

830 rain rate thresholds. Cases shown include those from Figure 3 along with intermediate 

83 1 rain thresholds. 
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TES/MLS CO at 562 & 215 hPa, TES akcl >= 0.080 










